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’’Orlando TDWR Testbed and Airborne Wind Shear Data Comparison Results” 
Steven D. Campbell, Anthony J. Berke and Michael P. Matthews 

Lincoln Laboratory 

Massachusetts Institute of Technology 
Lexington, Massachusetts, U.S.A. 

Telephone: 617/981-3386, E-mail: sdc@ll.mit.edu 


The focus of this talk is on comparing Terminal Doppler Weather Radar 
(TDWR) and airborne wind shear data in computing a microburst hazard index 
called the F factor. The TDWR is a ground-based system for detecting wind shear 
hazards to aviation in the terminal area. The Federal Aviation Administration will 
begin deploying TDWR units near 45 airports in late 1992. As part of this develop- 
ment effort, M.I.T. Lincoln laboratory operates under F.A.A. support a TDWR 
testbed radar in Orlando, FL. 

During the past two years, a series of flight tests has been conducted with 
instrumented aircraft penetrating microburst events while under testbed radar sur- 
veillance. These tests were carried out with a Cessna Citation II aircraft operated by 
the University of North Dakota (UND) Center for Aerospace Sciences in 1990, and a 
Boeing 737 operated by NASA Langley Research Center in 1991. A large data base 
of approximately 60 instrumented microburst penetrations has been obtained from 
these flights. 

The test flights in 1990 included the first-ever demonstration of real-time 
transmission of TDWR microburst graphical warnings to an aircraft for cockpit dis- 
play. A similar demonstration was carried out in 1991, with the TDWR microburst 
alerts being used to direct the NASA aircraft in making microburst penetrations. 

Post-flight analysis was performed under NASA funding to compare the F 
factor (Bowles & Targ, 1988) as measured by aircraft in situ sensors and estimated 
from TDWR microburst alarms. It was found that improvements are needed in the 


The work described here was performed under Air Force Contract No. F I9928-90-C-0002, 
and was sponsored by the Federal Aviation Administration and the National Aeronautics and 
Space Administration. The IJnted Stales ( u >vernment assumes no liability for its content or 
use thereof. 
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TDWR microburst alarm generation process to allow the aircraft F factor to be esti- 
mated accurately. These improvements include: shear-based outflow detection, 
physical model-based alarm representation, and compensation for the dependence 
of outflow intensity on altitude. The rationale for these improvements will now be 
discussed. 

The aircraft F factor can be estimated from TDWR microburst alarms using a 
formula proposed by Bowles (1988): 

F tdwr = K’ (AV/AR) | GS/g + 2h/TAS ] = F x + F z ( 1 ) 

where AV/AR is the TDWR-measured shear, GS is the aircraft ground speed, g is 
gravitational acceleration, h is the radar beam height and TAS is the aircraft true 
airspeed. K’ is a factor which attempts to relate the average shear in the microburst, 
AV/AR, to the peak shear in the microburst over a 1 km distance. The GS/g term 
corresponds to the horizontal component of F (F x ) and the 2h/TAS term is an esti- 
mate of the vertical (downdraft) component of F (F z ). It should be noted that the 
equation assumes that the aircraft penetrates through the center of the microburst. 

It was found that applying Equation 1 to current TDWR microburst alarms 
often overestimates the aircraft F factor. Examination of TDWR radar data shows 
that strong microbursts often contain small regions of intense shear inside a larger 
region of less intense shear. These intense shear regions are not identified by the 
current microburst detection algorithm, which attempts to identify the peak-to-peak 
velocity loss, rather than shear. Because of this, the shear associated with a micro- 
burst alarm is underestimated for strong microbursts. Applying the K’ factor to this 
underestimated shear leads to the correct F factor estimate for strong microbursts, 
but overestimates the F factor for weak microbursts. 

In order to better quantify the shear for use in Equation 1, a least-squares 
shear estimator was developed. The base polar radar data was first smoothed using a 
0.5 km x 0.5 km median filter. The least squares estimator was then applied over a 
seven-gate window ol TDWR velocity data for an effective distance of 0.9 km (i.e., 6 
gates eenter-io-center x 150 m per gate). The corresponding shear values were then 
applied to the following equation: 

F shrak - (dV/d R)| h [ GS/g + 2h/TAS | ( 2 ) 

where (dV/dR)| h is the least-squares shear at the radar beam height. 

ft was found that Equation 2 was an improvement but still often overestimated 
the aircraft F factor. Further examination of the radar data showed that there was a 
strong dependence of the outflow strength on altitude. Work by Mark Isaminger and 
Paul Biron of Lincoln showed that the outflow strength decreases linearly with height 
above the surface. This result was consistent with an analytical model of microburst 
outflows developed by Vieroy of NASA Langley (1991); this model is a modification 



of an earlier model developed by Oseguera and Bowles (1988). In the Vicroy and 
Oseguera & Bowles models, the horizontal shear is described by a shaping function, 
p(z), which is zero at the surface, reaches a peak at height h m and then drops off with 

increasing altitude. 

Using the altitude shaping function, p(z), the horizontal shear at the aircraft 
altitude, a, can be estimated: 

(dV/dR)| a = (dV/dR)|h I p(a)/p(h) | ( 3 ) 

and the revised F factor estimate can be written as: 

F alt.corr. - (dV/dR)|, ( GS/g + 2o/TAS | (4) 

where we now use the aircraft altitude, a, in (he downdraft estimation term, 2a/TAS. 
This formula reflects the concept that as the aircraft altitude increases, the horizontal 
shear will decrease but the downdraft component will increase. 

Equation 4 was found to estimate the F x component quite accurately, but still 
tends to overestimate the F z component. Further reflection shows that the 2a/TAS 
term leads to an overestimate of the vertical component, since it is assumed that t e 
aircraft flies directly through the center of the microburst. In fact, many of the pene- 
trations were made at the edge of the outflow where the Vicroy model predicts an 
updraft, rather than a downdraft. 

Accordingly, a final modification was tested which divided the aircraft data 
into center and edge penetrations. For center penetrations, the unmodified Equation 
4 was used; for edge penetrations, the vertical component estimator was changed 
to -a/TAS (i.e., an updraft at the edge equal to half the center downdraft): 

F hor.corr. = (dV/dR)la | Ofi/g ♦ 2a/TAS I, center ( 5a ) 

= (dV/dR)j ; , | GS/g - a/TAS |, edge ( 5b ) 

Applying Equation 5 yielded an improvement in the mean F z component, however, 
the data points were clustered as either loo high or too low. A further refinement 
would be to scale the vertical compensation according to distance Irom the outflow 

center. 

These results lead to the notion that several improvements could be made to 
the existing TDWR microburst recognition algorithm to allow accurate F factor esti- 
mation. First, shear-based outflow detection at multiple thresholds would allow re- 
gions of intense shear to be identified inside of larger outflow regions. Second these 
shear regions could be used to create a microburst representation based on a physica 
model consisting of an outflow center and an outflow edge. Third, an analytic micro- 
burst model or other technique could be used to compensate for the dependence o 
outflow intensity on altitude. Fourth, the improved microburst representation could 
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be used to estimate the vertical component of the microburst based on distance from 
the outflow center. 

A key goal for operations during the summer of 1992 will be to more accurate- 
ly characterize the altitude dependence of microburst outflows. It is planned to ac- 
complish this goal by carrying out rapid, low-altitude scans of microburst outflows 
by three radars during aircraft penetrations. The three radars will be the TDWR 
testbed plus two C-band radars operated under F.A.A. funding by the University of 
North Dakota and Massachusetts Institute of Technology. These radars are situated 
in such a fashion to allow triple-Doppler reconstruction of the three-dimensional 
wind fields at the Orlando airport. These triple-Doppler wind field reconstructions 
will allow both the horizontal and vertical components measured by airborne and 
ground-based sensors to be compared. 

In summary, a large data base of instrumented microburst penetrations while 
under TDWR testbed radar surveillance has been obtained over the past two years at 
Orlando. These tests also marked the first-ever demonstration of real-time data link 
transmission of TDWR microburst alerts to aircraft for graphical display in the cock- 
pit. Additional flight tests will be performed in 1992, including penetrations with 
rapid update, low-altitude triple-Doppler radar scans. 

Sixty microburst penetrations have been examined to determine how well the 
aircraft F factor can be estimated from TDWR data. Analysis of the data shows that 
several improvements to the current microburst recognition algorithm would be 
needed to allow the aircraft F factor to be accurately estimated. These improvements 
would improve the quality of the microburst alerts currently supplied to ATC person- 
nel and, in the future, supplied to pilots directly via Mode S Data Link. 

References: 


1. Bowles, R.L. and R. Targ: Wind Shear Detection and Avoidance: Airborne 

Systems Perspective, International Congress of Aeronautical Sciences, 
Jerusalem, Israel, August-September, 1 088. 

2. Vicroy, D.D.: A Simple, Analytical, Axisymmctric Microburst Model for 

Downdraft Estimation, NASA Technical Memorandum 104053, NASA Langley 
Research Center, Hampton, VA, February, 1991. 

3. Oseguera, R.M., and R.L. Bowles: A Simple, Analytic 3-Dimensional Downburst 

Model Based on Boundary Layer Stagnation Flow, NASA Technical Memo- 
randum 100632, NASA Langley Research Center, July, 1988. 


816 



OUTLINE 



• SUMMARY 




TERMINAL DOPPLER WEATHER R 


< 

DC 

O 

o 

DC 

Q. 

DC 

< 

O 

< 


CM 

O 



O 

O) 


CD 

CO 


CO 

CO 


r- 

oo 


CO 

CO 


in 

co 


if) 

h- 

Z 

LU 


_ DC 

Q LLJ 

a. 

LU 5 

Li uj 



cd 


o a 
r z 

* < q 

0 CD ^ 


D o 

o 

LU 

CD 

H 

if) 

LU 


< 

CD 

CD 

DC 

X 


o 

UJ 


z 

UJ 


CD 

> 



(D 

z 


Z 

o 

UJ 

h- 

(J) 

5 

DC 

3 

CD 

o 

DC 

H 

O 

5 

CD 

< 

u 

a 

< 

X 

3 

^CD 
Q > 

X 

n 

u. 



z 

o 

< CD 

i- 

X 

DC 

a 

h 

CD 

3 

2 

X 

o 

.j 

X 

O 

h- 

o 

X 

cc 2 

X 3 
UJ 

o 

5 

o 

H- 

X 


X 

z 

o 



CD 

$ 

Q 


< 






H 


< 


< 






u. 




s 


818 


TDWR (45 SYS) 
RAYTHEON 





WINDSHEAR DETECTION II 
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ORLANDO ’91 TDWR TESTBED ACTIVITIES 
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WIND SHEAR WARNING PROVIDED OPERATIONALLY 
TO ORLANDO ATC DURING 6 WEEK PERIOD 




FLIGHT TEST ACTIVITY AT ORLANDO 
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DATA LINK TO NASA B737 
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NASA COCKPIT DISPLAY 




TDWR F FACTOR ESTIMATE 
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CONTINUITY 
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MICROBURST REPRESENTATION 
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WEAKER SHEAR & UPDRAFT AT EDGE 
ALTITUDE DEPENDENT 




MICROBURST REPRESENTATION 
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- CENTER NOT WELL LOCALIZED 

- EDGE EXTENT UNDERESTIMATED 

- NO ALTITUDE DEPENDENCE 



MICROBURST REPRESENTATION 
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- SHEAR-BASED OUTFLOW DETECTION 

- LOCALIZE CENTER AND EDGE 

- ALTITUDE COMPENSATION 


SHEAR COMPUTATION 
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Mean Vertical Velocity Structure of Orlando Microburst 

Delta V Max = 22m /$ 



DeltaV (m/s) 






Vertical Shear Structure of Orlando Microburst 
6/20/91 21:24 



Shear (m/s/km) 
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HORIZONTAL SHEAR VS. ALTITUDE 
(VICROY/O&B MODEL) 




TDWR Shearmap vs. Aircraft Total F-Factor 
using Vicroy Model Correction 



Plotting Time => 04/07/1992 18:16:16 








TDWR Shearmap vs. Aircraft Horizontal F-Factor 
@ Ftotal Peak Time using Vicroy Model Correction 


o 



.lojorj-j pojBiupsg yA\a.I, 


839 


Plotting Time => 04/09/1992 14:17:29 




TDWR Shearmap vs. Aircraft Vertical F-Factor 
@ Ftotal Peak Time using Vicroy Model Correction 
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Plotting Time => 04/10/1992 11:37:51 



HORIZONTAL OFFSET COMPENSATION 
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- POSITIVE 2a/TAS FOR CENTER 

- NEGATIVE a/TAS FOR EDGE 

- NEED TO LOCALIZE OUTFLOW CENTER & EDGE 
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OPERATIONAL TESTS IN 1993-1994, FOLLOWED BY 
INITIAL OPERATIONAL CAPABILITY (IOC) IN 1996 




SUMMARY 



845 


RAPID TRIPLE-DOPPLER LOW-ALTITUDE SCANS 



The Orlando TDWR Testbed and Airborne Wind Shear 
Data Comparison Results 

Questions and Answers 

Q: Dan Vicroy (NASA Langley) - You pointed out some improvements or possible 
improvements to the TDWR algorithms. Can you comment on the implementation issues and 
what kind of time line you are looking at for implementing these improvements? 

A: Steve Campbell (MIT Lincoln Lab.) - The TDWR was implemented as a very fast track 
program. We knew that there would be some refinements. When the TDWR was designed, the 
idea was that all you needed to do was detect the change in velocity. I think we now understand 
that it is not true. There are really two avenues through which we could make improvements. 
One is that the FAA expects to upgrade the TDWR algorithms over a period of time. The other 
is that there is another program which is starting up called the Integrated Terminal Weather 
System Program in which we will be incorporating data from a number of sources, TDWR, 
surface observations and aircraft data. That may also be an avenue for making these 
improvements. As far as how long that is going to take, well it is going to take some years. I 
think we are at least plugged into that process. 


846 



